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We study a microscopic mean-field model for the dynamics of the electron glass, near a local 
equilibrium state. Phonon-induced tunneling processes are responsible for generating transitions 
between localized electronic sites, which eventually lead to the thermalization of the system. We 
find that the decay of an excited state to a locally stable state is far from being exponential in 
time, and does not have a characteristic time scale. Working in a mean-field approximation, we 
write rate equations for the average occupation numbers (jit), and describe the return to the locally 
stable state using the eigenvalues of a rate matrix, A, describing the linearized time-evolution of 
the occupation numbers. Analyzing the probability distribution P{\) of the eigenvalues of A we 
find that, under certain physically reasonable assumptions, it takes the form P(A) ~ jjj, leading 
naturally to a logarithmic decay in time. While our derivation of the matrix A is specific for the 
chosen model, we expect that other glassy systems, with different microscopic characteristics, will 
be described by random rate matrices belonging to the same universality class of A. Namely, the 
rate matrix has elements with a very broad distribution, i.e., exponentials of a variable with nearly 
uniform distribution. 
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PACS numbers: 71.23.Cq, 73.50.-h, 72.20.Ec 



I. INTRODUCTION 

Experiments conducted on thin-films of amorphous or 
crystalline semiconductors such as indium-oxide or sili- 
con, show that when driven out-of-equilibrium (for ex- 
ample by shining light on the system or by changing a 
gate voltage), the system exhibits slow relaxations, ob- 
servable on the scale of minutes or hours^. In many 
cases a logarithmic or weak power-law time-dependence 
of the measured quantity (such as conductance and ca- 
pacitance) is observed over many decades of time^£. A 
common feature of the experimental systems is that they 
are highly disordered, so that most electronic states are 
localized. If the carrier concentration is high enough 6 , 
the (unscreened) Coulomb interactions may play an im- 
portant rolei. This system is usually referred to as the 
electron glass, since it exhibits many features characteris- 
tic of glassy systems: memory effects— (the conductance 
depends on the previous perturbations applied to the sys- 
tem) and aging^ (the duration of time of the perturbation 
is applied affects the relaxation timescales). Similar ef- 
fects have been observed in granular A l 10 ' 11 showing that 
the underlying principles may be more general. 

In this paper we study a mean-field model for the dy- 
namics of the system. A variety of systems in nature 
can be described, near a locally stable state, by a matrix 
equation of the type: 



dSn A -> 
— — =A-5n, 
dt 



(1) 



the component Srii = n, — /; is the deviation of the av- 
erage occupation of the i'th site, n*, from its value ft 
at the locally stable point. The local stability of the 
point implies that the matrix A must have only non- 
positive eigenvalues, and, for large systems, their distri- 



bution will determine the average time-dependence of the 
return to the locally stable point, after the system was 
slightly pushed away from it. It must be emphasized that 
our approach is different from the usual theoretical ex- 
planations of aging phenomena in glasses, in which the 
system explores the energy landscape, and slow relax- 
ations are a result of the existence of many mctastablc 
states. In our model, the system is found in the vicinity 
of one locally stable point, at all times (we do not use the 
term metastable to stress this difference). This assumes 
that the initial perturbation is small enough (and so is 
the temperature), such the system does not reach other 
(lower) minima, but remains in the same region of phase 
space. Slow relaxations are due to isolated states that, 
statistically, happen to have a long life time. It should 
be emphasized that although the interactions lead to the 
non-trivial Coulomb gapi£ in the equilibrium state, the 
slow dynamics will occur also without interactions. 

If it is given that the distribution of eigenvalues di- 
verges at small (negative) eigenvalues, and is of the form 
P(X) ~ ( as happens in our model) it is straightfor- 
ward to see that logarithmic relaxation in time, in an ap- 
propriate time window, is obtained (assuming the eigen- 
vectors are excited with uniform probability). In this 
work we show that starting from a realistic microscopic 
model for the electron glass system, the described situa- 
tion indeed occurs, and we argue that it is plausible that 
other physical systems will also show similar results. 

The structure of the paper is as follows. The model 
is defined in III Al In III B( we review the application of 
the mean-field approximation to the peculiar equilibrium 
properties of the system, manifesting the Coulomb gap. 
In III CI we briefly discuss the mean- field steady-state so- 
lution in the presence of an external field, leading to the 
Miller- Abrahams models. 



2 



In a similar fashion, in section \UJ\ wc suggest to study 
the dynamics of the system by writing a set of ordinary 
differential equations, described by Eqs. |J2J) and ((6]), giv- 
ing the time-evolution of the occupancies of the localized 
states. This is already an approximation neglecting in- 
terference or quantum fluctuation effects. In IIII A[ we 
study the dynamics of the electron glass, starting from 
an out-of-cquilibrium state. Linearizing Eq. we ob- 
tain the time-evolution equations of the occupations and 
obtain Eq. ([T]), with the random matrix A belonging to 
a different class from the usual gaussian random matrix 
ensembles. The statistics of the eigenvalues is studied 
numerically, see Fig. [3] In IIII Bl wc study a simplifying 
limit, analytically. Both lead to a distribution of eigen- 
values P(A) diverging at low values (down to a cutoff), 
leading to slow relaxations of the physical observables, as 
seen experimentally. This behaviour might be character- 
istic of glassy systems. Finally, in IIII CI we discuss the 
relation between the relaxation of the occupation num- 
bers and the conductance. 



II. MEAN-FIELD MODEL FOR ELECTRON 
GLASS 

In this section we discuss a specific microscopic model 
for the dynamics of the electron glass. We will show that 
it leads to a rate equation of the type of Eq. |T]), and 
find explicitly the matrix elements. 



7ij ~ |M,| a i//,(l - fj)e—t [1 +N(AE)}, (2) 

where /{ is the Fermi-Dirac distribution. For upward 
transitions (Ej > Ei) the square brackets are replaced 
by N(AE). These rates may be rcnormalized due to 
polaron-type orthogonality effects^. 

We will be interested in the dynamics of the sys- 
tem when it is out-of-equilibrium, namely in the time- 
dependence of the occupation numbers and the conduc- 
tance after an initial excitation. But we first show how 
non-trivial equilibrium properties are obtained from the 
mean-field picture. 



B. Equilibrium properties near a locally stable 
point 

In an approximation similar to those used in spin glass 
theory^, we define /, = (rii), where n, the site occupa- 
tion (which takes the values zero or one) and () denotes 
averaging over a time scale much larger than that related 
to the frequency of the phonon processes but smaller than 
the relevant scale for the observation of the dynamics. 
This is a mean-field type approximation, and may be 
used regardless of the interactions in the system. Let us 
first discuss the thermal equilibrium state, near the lo- 
cally stable point. The sites occupation must follow the 
Fermi-Dirac distribution and therefore: 



A. Definition of the model 

We study a system of N localized states and M < N 
electrons, with a coupling between the electrons and a 
phonon reservoir. Since the states are localized, the elec- 
trons will interact via an unscreened Coulomb potential. 
In the absence of electron-electron interactions, the lo- 
calized states have different energies, ej, due to the dis- 
order. Our model also contains structural disorder: the 
positions of the sites are assumed to be random. Al- 
though localized states are orthogonal, their tails over- 
lap, and therefore phonons may induce transitions be- 
tween them. The generic coupling between electrons and 
phonons is given by the form ^ M q c\cj (b q + b- q ) where 

cj , Cj are electron creation and annihilation operators at 
local sites i,j and b q annihilates a phonon. M q is a coef- 
ficient accounting for the strength of the electron-phonon 
coupling. 

Let us denote the energy difference of the electronic 
system before and after the tunneling event by AE, con- 
taining the interaction effects. For weak electron-phonon 
coupling (\M q \ 2 v -C AE, where v is the phonon density 
of states), the transition rate 7y of an electron from site 
i with energy Ei to site j with energy Ej < Ei a distance 
Tij away, can be calculated treating the coupling as a per- 
turbation. This yields, up to polynomial corrections-^: 



fi(Ei) = 



1 + e" 



(3) 



where fi is the chemical potential, and the Boltzmann 
constant is set to be one. 

In the mean-field approximation we can calculate the 
average potential energy of site i: 



E, 



... If) 



(4) 



This approximation improves as the number of inter- 
acting sites increases^. The long range nature of the 
interaction means that the energy of a site will be de- 
termined by many of its neighbours, and gives intuitive 
justification of the use of mean-field theory. Neverthe- 
less, we should emphasize that this is an uncontrolled 
approximation, and the limits of its validity should be 
checked. 

Combining Eq. ([3]) and one obtains a self- 

consistent equation for the energies. It is common to use 
an unbiased disorder distribution, and add a background 
charge v to each sit©i£. In the mean-field picture, this 
will lead to the equation: 



Ei 



E 



e t 



(5) 
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For half filling, (j, = 0, 



0.5. 



One should notice that there are many solutions to 
this self-consistent equation. Rigorously, one cannot call 
any solution an equilibrium distribution, since the equi- 
librium distribution is a Boltzmann average over all con- 
figurations, not only those near the locally stable point. 
The solution may be viewed as a 'local equilibrium'. We 
will see that the physical picture obtained is quite plausi- 
ble. At low temperatures^, the probability distribution 
of the energies will contain a soft gap at the Fermi energy, 
known as the Coulomb g a plilL12JlUL2i212i. 



Eq. (|5|) can be solved numerically by starting with 
a random set of energies, and evolving them iteratively, 
within the mean-field model. This was done following 
Ref. [l^ . by solving the equations for many random 
instances, and averaging over them. In this way a his- 
togram of the on-site energies is obtained. When nor- 
malized correctly, it gives the single-particle DOS (den- 
sity of states), as function of energy. The results for 
two-dimensions, yielding the Coulomb gap, are given in 
Fig. [TJ Notice that the obtained DOS contains a lin- 
ear gap near the Fermi energy, in accordance with other 
work s 19 ' 23 . 




FIG. 1: Histogram of the site energies in two dimensions, for 
N=10000, fx = (half filling), obtained by solving the self- 
consistent Eq. Ipa). The sites were uniformly distributed in 
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a square, with - e T = 20, where r nn denotes the average 
nearest-neighbor distance and T is the temperature. The en- 
ergies £i were uniformly distributed in the interval [— 
with W = 1. The y axis denotes the probability density of 
the energies Ei. The graph is the average over 300 instances. 
Notice the finite value of the density at the minimum, due to 
the finite temperature. 



Response to an external field: steady-state 
solution 



When a small electric field is applied, there are cor- 
rections to the average occupations and also to the av- 
erage energies. It can be shown that the problem of 
finding the steady-state solution corresponds to that of 
solving the steady-state of a resistance network, using 
Kirchoff's laws^ 3 -. The solution, when neglecting in- 
teractions, gives the well-known Mott Variable Range 
Hopping 2 ^, which was observed experimentally in many 
cases^. We emphasize that this calculation is in fact a 
mean-field one: the steady-state solutions are obtained 
from timc-dcpcndcnt equations which arc essentially the 
mean-field equations. In the following we propose to use 
the same ideas to discuss the dynamics of the system 
out-of-equilibrium. 



III. THE DYNAMICS 

Let us pose the following question: How will the oc- 
cupation numbers or conductance depend on time, when 
the system is pushed slightly out of the locally stable 
point? Having seen that the mean-field approximation 
yields the correct density of states as well as the out-of- 
cquilibrium steady-state solution in the presence of an 
external field, we propose to use the same approximation 
to describe the dynamics of the systems, when prepared 
out of local equilibrium. 

Experimentally, the form of the relaxation depends on 
the details and mechanism of the excitation. For simplic- 
ity, let us assume that the initial perturbation takes the 
form of a random addition drii to the state occupations, 
with $ n i = 0, reflecting particle number conservation. 

Assuming that the initial change in the occupations is 
small, wc can still use Eq. @ for the tunneling rates, 
with the average occupations at the locally stable point 
fj substituted by the occupation numbers slightly out of 
equilibrium rij (which can take any value between and 
1). The energies at each instance are related to the out- 
of-equilibrium occupations by Eq. ^ , upon replacing fj 
by rij , and we can write the time-evolution of the average 
occupation as: 



(6) 



This defines the problem completely. At the locally 
stable point itself, the RHS of Eq. ([6]) must vanish. 
Therefore not too far from the point, we can take the 
first (linear) order in the quantities 5rii, the deviations 
from the stable point. The linearized equation then takes 
the form of Eq. ([1]) , where 8n is a vector of the deviations 
of the occupation numbers from their local equilibrium 
values. The N eigenvalues and eigenvectors of the matrix 
A will determine the decay rates of the system. 
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We would like to stress that the matrix A can be calcu- 
lated for the cases of interest, by linearizing the equations 
of motion near the locally stable point. This strategy 
is completely general, and will be valid for any system 
which can be described by equations of motion, and has 
a locally stable point (such would be the case for most 
classical systems, and many quantum systems in a mean- 
field approximation). The dynamics of the system when 
pushed slightly away from the fixed point, will be charac- 
terized by the eigenvalues of the rate matrix. In the com- 
mon case where disorder plays a role, the dynamics will 
depend on the distribution of eigenvalues of the matrix, 
thus, we obtain a problem of random matrix theory^, 
where the eigenvalue distribution is responsible for the 
dynamics of the system. An extremely relevant prop- 
erty of the electron glass case, as we shall demonstrate 
in IIII A[ is that the entries of the random matrix are ex- 
ponentials of the broadly distributed parameters (energy 
and distance) . Another important feature of these matri- 
ces is that the sum of every column vanishes. These prop- 
erties make this matrix belong to a different class from 
the usual gaussian ensembles treated in random matrix 
theory, and will play an important role in the dynamics, 
leading to slow relaxations. A similar class of matrices 
was studied previously by Mezard et al.— . 

In the following section we will derive the form of the 
matrix A for the particular case of localized states cou- 
pled due to phonons. We will find that the probabil- 
ity distribution is divergent for small eigenvalues, and 
suggest what the minimal properties leading to such a 
distribution are. The implications of this distribution on 
the time-dependent relaxation of the occupation numbers 
and conductance will then be discussed. 



A. Application to the electron glass model 

Starting from Eq. © , a calculation of the elements of 
matrix A in Eq. © shows that: 



exponentially, meaning the 



term explodes much 



7 P. 



(i- nj y 



(7) 



where 7? are the local equilibrium rates, given by Eq. 
©. For i^j: 



1 



Ujil — nj) 



(8) 



Notice that the matrix is not symmetric, due to the 
Uj(l — nj) term. The sum of each column of the matrix 
vanishes, guaranteeing particle number conservation. 

A-priori one would expect that at low enough tem- 

2 

peratures T <C , we could neglect the first term in 
the equation for the regime of interest. However, at low 
temperatures the occupations of the sites tend to or 1 



faster than the i part in the second term 



Uj (1— Uj) 

Viewed in a 

different way, if one looks at the expression of the mean- 
field rates 7^ ~ n^(l — n,j)[l + N(AE)], one sees that if 
two states are close in energy, then the first term in the 
matrix element flS) ~ ^-N(AE) ~ -^ E . Therefore there 
is good coupling between any two states close in energy 
(and distance), not only those ones close to the Fermi 
level, as is the case for the second term. Therefore the 
'phase space' is much larger for the first term, and the 
second one can be neglected. We have calculated numer- 
ically the eigenvalue distribution for some specific system 
parameters, and indeed it was found that the second part 
has a small influence on most eigenvalues, see Fig. [2j 

An important property follows: the off-diagonal ele- 
ments are positive. Together with the property that the 
sum of every column vanishes, the stability of the mean- 
field solution is guaranteed: all the eigenvalues are nega- 
tive, characterizing decay. A proof of the statement can 
be found in Ref. [28l | . 

Let us consider the distribution of the eigenvalues. Fig. 
Oshows the distribution of eigenvalues of the matrix A, as 
obtained numerically. We first found a mean-field solu- 
tion by iterating the equations (see, for example^), then 
used Eq. |8]) to construct the relevant matrix. The eigen- 
values of this matrix were then found numerically Notice 
that the localization length £ influences the dynamics, al- 
though it has no effect on the equilibrium properties, at 
least as long as the localized states are spatially well scp- 
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FIG. 2: Comparison of the eigenvalues of the full linearized 
matrix (stars) and the ones obtained after neglecting the 
'Coulomb' term (crosses), the second one in Eq. JS). The 
difference is mostly seen for the large magnitude eigenvalues. 

= 10. The energies ti were 



TV = 1000, = 0.1 and 



uniformly distributed in the interval M-], with M = 10. 

The inset shows the relative error in replacing the full matrix 
with the approximated one, defined as the difference between 
the approximation and exact diagonalization, divided by the 
exact value. 
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arated. 

In IIII Bl wc will analyze a simplifying limit, when the 
rather complicated dependence on energy can be ne- 
glected, and consider only the exponential dependence 
of the tunneling rate on length. Both limits give approx- 
imately a -pjy distribution (up to logarithmic corrections), 

reminiscent of j noised. This suggests that the result 
may be more general, and not dependent on the details of 
the specific model. Note that the interactions affect the 
mean- field solution (and the Coulomb gap), but the cal- 
culation shows the slow dynamics will exist also without 
them. 

Let us now discuss the consequences of this distribu- 
tion for the dynamics of the system. Having found that 
the distribution is ~ jL down to some minimal value, 

we conclude that if all eigenvectors (except the one with 
eigenvalue zero, not conserving particle number), are ex- 
cited with equal probability, the time-evolution of the de- 
viation from the local stable point will be related to the 
Laplace transform of the eigenvalue distribution. This 
will give rise to logarithmic decays. Let us show this in 
more detail: if the eigenvectors are denoted by ri\, the 
time-dependent deviation from the locally stable point 
will be given by: 



5n = J2 c ^ e ~ M - ( 9 ) 

A 

We shall assume the eigenvectors are excited with 
roughly uniform probability. Going to the continuous 
limit, and utilizing the fact that the components of the 
eigenvectors themselves are also random variables, we ob- 
tain that the norm of Sn should relax as: 




2 4 6 8 

X x10" 6 

FIG. 3: Distribution of eigenvalues obtained numerically, av- 
eraged over 1000 realizations. The parameters were as for 
Fig. [3 The fit to the reciprocal distribution is linear. 



\6n(t)\~ -IE- \og[t\ mm ] , (10) 

for v^— < t < T J— . 



B. Dynamics of exponential models 

Hitherto we have discussed a specific model for glass 
dynamics in a system constructed from interacting elec- 
trons and phonons. The actual form of the rate matrix 
eigenvalue distribution did not depend strongly on the 
details of the matrix elements. In this section we will 
show that there are few sufficient conditions on the ran- 
dom rate matrix A that will make the relaxation process 
long. 

Let us look at the dynamics which follow from a class 
of random matrices obeying the following properties: 

1. The sum of every column vanishes. This follows 
from particle number conservation. 

2. The entries of the matrices arc distributed over a 
very broad range. This happens, for example, when they 
are exponentials of a more or less flat distribution^. 

We expect that a variety of systems that exhibit a 
glassy behaviour may be described by a random rate ma- 
trix belonging to this class. The matrix obtained for the 
electron glass system indeed obeys these properties. The 
first property was shown explicitly in section IIII Al To 
see the second, let us examine Eq. (J8j) - If we neglect cases 
where the energies Ei, Ej and \Ei — Ej \ are smaller than 
T, we can recast the equation into a more transparent 
form: 

A i:j ~e f e . (ii) 

Due to the exponential, the matrix entries are indeed 
broadly distributed. 

We shall now discuss a specific class of matrices which 
can be analyzed analytically. As seen in Eq. (fTTj) , the 
matrix elements for the electron glass system contain a 

factor e « . If £ is much smaller than the typical dis- 
tance r„„, it is plausible that this factor would be dom- 
inant in determining the eigenvalue distribution. This 
motivates us to discuss a simpler model, of so-called 
distance matrices^: assume we have TV random points 
in a two dimensional space. Let us define a matrix 

Bij = e « , where r,j is the distance between points 
i and j, and £ some constant. Let us choose the diagonal 
elements of the matrix such that the sum of every col- 
umn vanishes. Following the previous discussion of the 
dynamics, we are interested in the distribution of eigen- 
values of such a matrix. A mapping of this problem to 
a field theory problem is given in Ref . [28[ , enabling one 
to look at a low-density approximation to the theory. 
Mezard et. al calculate the resolvent R = jrTr . , the 
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imaginary part of which yields the density of states, i.e, 
the distribution of eigenvalues. Using their formula (21) 
for the case of f(r) = e~ « , we obtain that the low-density 
expansion of the resolvent is: 



the probability that its nearest-neighbour is at least dis- 
tance r away. This is equivalent to asking that all of its 
neighbours are at least a distance r away, and since they 
are randomly distributed, we obtain that: 



R{\) 



wj dxdy { 



1 



X + 2e~ 



(12) 



where the integrals are performed over the whole volume. 
The second term gives rise to a delta function at the 
origin, which comes from the zero eigenvalue the matrix 
always possesses, and is of no particular interest since 
the eigenvector associated with this eigenvalue cannot 
be excited while preserving the particle number. The 
condition for the approximation to be valid is £ -C r nn , 
as we shall show later in a more transparent way. 



Since the density of states is given by — 



Im[R(,X)] 



we 



can use the fact that Im\—r-A 
DOS as: 



-inS(x) and obtain the 



P(X) = 
V ; 2V 



dxdyd(X + 2e «)• 



(13) 



Performing the integral in one-dimension, for eigcnval- 
ues not too close to the minimal value 2e « , leads to the 
result: 



P(X) = 



LX ' 



(14) 



with A in the interval [—2, — 2e~]. 

Repeating the calculation in two-dimensions, again, for 
eigenvalues not too close to the minimal values, yields: 



P(A) 



irN£ 2 log(- 



L 2 X 



(15) 



We shall now give a transparent demonstration of these 
results. In the low density limit, we can couple each site 
to its nearest-neighbor, thus dividing the system into ^■ 
pairs. Neglecting the effect of all other sites, we obtain 
that the eigenvalues will be similar to those of an ensem- 
ble of 2X2 matrices of the form: 



M 




(16) 



Since the matrix elements of B decay on the scale of 
£, this would be a good approximation for £ <C r nn . One 

— I a- — y\ 

eigenvalue of M is 0, and the other is — 2e « . There- 
fore half the eigenvalues of B will be the zero under this 
approximation, and the distribution of the other eigen- 
values will be that of the random variable — 2e~ , where r 
is the nearest-neighbour distance. Notice the zero eigen- 
values correspond to the second term in Eq. (|12D . 

The distribution of the nearest-neighbour distance can 
be calculated: looking at a typical site, let us calculate 



Prob{r) = {l-V DTn ) N -\ 

where V\ = 2 and V% = ir. 

For r « I we can approximate: 



(17) 



Prob(r) = e 



(18) 



We have assumed that the initial site is not too close 
to the boundaries. Since we are interested in the proba- 
bility, the sites near the boundary will give a negligible 
correction to the above probability: the sites for which 
Eq. (Ti~8"|) fails are a distance of order r nn or less from 
the boundary. Therefore their fraction in the system is 
of order ^-j^. For N ^> 1, this fraction is negligible. 

The probability distribution can be calculated by dif- 
ferentiating with respect to r, leading to: 



P{r) = V d DN- 



,D-1 



L L 



(19) 



By construction, the probability distribution is exactly 
normalized. 

In one-dimension, the eigenvalue distribution that fol- 
lows is: 



e — L-^l'ogC— SOI 
P(A) = -N £ 



with e = ^y-S <c 1, while for two-dimensions: 



(20) 



P(A) = 



irNe log(-|)e- 



r Wlog 2 (-|) 



L 2 X 



(21) 



Aside from the exponential term, Eqs. (|20|) and (|2"Tj) 
coincide with the field-theory Eqs. (fT4|) and (fT5|) . No- 
tice that in the latter there is a cutoff on the eigenvalue 
magnitude, while for Eqs. (Ti~4")) and (|15p the distribu- 
tion is nonzero also for very small eigenvalues. Fig. [4] 
shows the results of numerical simulations for the case of 
low-density distance matrices in two-dimensions, and a 
comparison to the theory. 



C. Relaxation of the conductance 

In many cases a logarithmic relaxation in time is ob- 
served for the conductanc e) 3 ' 4 ' 11 . One should ask what is 
the relation between the relaxation of the conductivity 
and that of the occupations, for the electron glass model. 
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We now present an intuitive argument motivating the 
speculation that the time-relaxation of the conductance 
should be similar to that of the occupations. The essence 
of the argument is the claim that any perturbation of the 
equilibrium configuration will lead to enhanced conduc- 
tance: if this is true, it is reasonable that as the typi- 
cal deviation of the occupation number relaxes, so does 
the enhanced conductivity, until it reaches its equilib- 
rium value. For small enough deviations, the two will be 
proportional to each other, as one can always take the 
lowest order term in the expansion of the dependence 
of the out-of-equilibrium conductivity on the occupation 
number deviation. 

Let us explain our claim that any perturbation will 
lead to enhanced conductivity. This may come about by 
two physically different mechanism: first of all, we note 
that when the system is excited, we create vacant sites 
well below the Fermi surface, and add electrons above it. 
Electrons will tunnel between these sites, and thus even 
at very low temperatures current may flow through the 
system. 

The second mechanism is more subtle, and is related to 
the Coulomb gap. Let us look at the Einstein relation^, 
a = e 2 ^D. We do not expect to have any anomalies in 

the thermodynamic DOS, jp, but the diffusion constant 
D should be much smaller. This is because the single- 
particle DOS at the chemical potential vanishes: moving 
an electron from a site with energy close to the chemi- 
cal potential to another site, will necessitate an energy 
of order of the width of the Coulomb gap. Therefore 
the Coulomb gap significantly lowers the conductivity. 
We shall assume that for a finite, large enough, single- 
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FIG. 4: Distribution of the eigenvalues in two dimensions, 
for a low-density = 0.05). N = 1000, and the results 

were obtained after averaging over 1000 instances. The points 
were distributed uniformly in the unit square, and the matrix 

element Aij = e e , with £ = 0.0016. The dashed curve is a 
plot of Eq. (|15[) . while the solid curve is a plot of Eq. (|21[) . 



particle DOS at [i, the conductivity increases with the 
single-particle DOS at (for systems close to a local 
equilibrium, which exhibit the Coulomb gap). The tem- 
perature should be low enough such that the local equi- 
librium Coulomb gap would not be smeared^. Let us 
suppose that due to our initial perturbation of the lo- 
cal equilibrium configuration, we have some excess (pos- 
itive or negative) 8rii in the occupation number at site i. 
Assuming these numbers to be random with a standard 
deviation 8n, the energy at site j will now have an addi- 
tional contribution A finite single-particle DOS 
at /i will arise, proportional to bn. Both mechanisms 
show that within this model the conductance relaxation 
should have a similar time-dependence to that of the oc- 
cupation number relaxation. 



IV. CONCLUSIONS 



We have studied a finite temperature mean-field model 
for the dynamics of the electron glass system. For a per- 
turbation which drives the system not too far from the 
(mean-field) local equilibrium, we mapped the problem 
onto rate equations with a random relaxation matrix A. 
The matrix A belongs to a class different than the gaus- 
sian random matrix ensembles. We found that the dis- 
tribution of the eigenvalues is approximately ~ pjy, and 
naturally yields a logarithmic decay of the occupation 
numbers. This may lead to a logarithmic decay of the 
conductance. Such a logarithmic decay of the conduc- 
tance is observed experimentally in many cases. We em- 
phasize the remark made before that the j distribution 
of decay eigenvalues should be much more general than 
for the specific model considered. It might also hold, for 
example, in the case of multi-particle transitions, which 
are believed to be relevant for the long time properties of 
glasses. Further research is needed to obtain additional 
predictions of this model, such as the time-dependence 
of the Coulomb gap, and the voltage-dependent conduc- 
tance in the 'two-dip' experiment^. 
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